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Abstract. - 

We show that the density and temperature dependences of the a-relaxation time of sev- 
eral glassforming polymers can be described through a single scaling variable X — e(p)/T, 
where e(p) is well fitted by a power law p x , x being a species-specific parameter. This implies 
that "fragility" is an intrinsic, density-independent property of a glassformer characterizing its 
super-Arrhenius slowing down of relaxations, and it leads us to propose a modification of the 
celebrated Angell plot. 



The glass transition of liquids and polymers is conventionnally studied at constant, usually 
atmospheric, pressure (P) by cooling the system. However, in an effort to disentangle the 
effects of density (p) and temperature (T) that both influence the viscous slowing down under 
isobaric conditions and to provide a more stringent test of existing models and theories, there 
has recently been an increasing number of sytematic studies of glass formation using both 
T and P as experimentally controlled variables. Glassforming polymers are interesting in 
this respect : first, because of their rather high compressibility allowing to cover a significant 
range of densities and second, because of the important role played in polymer science by the 
free- volume theory that puts the emphasis on congestion effects due to the increase of density 
as one approaches the glass transition [1] . 

In this letter we analyze new dielectric relaxation data obtained on several glassform- 
ing polymers up to 3 kbars poly(epichlorhydrinc) (PECH), poly(vinylmethylether) (PVME), 
poly(vinylacetate) (PVAc), poly(methylmethacrylate) (PMMA)), and a combination of di- 
electric, neutron scattering and calorimetric data on 1,4-poly (butadiene) (1,4-PB). We show 
that the p and T dependences of the a-relaxation time (r Q expressed, say, in seconds) can be 
described by a single scaling variable X = e(p)/T, 

log(T a (p,T))=F(^), (1) 

where, over the range of densities experimentally accessible, e(p) can be fitted by a power 
law, p x ; the parameter x and the function F(X) are species-specific. This scaling generalizes 
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to polymers work recently done on several molecular liquids. As a spin-off, we show that 
"fragility", a concept introduced by Angell [2] to characterize the degree of deviation from 
Arrhenius behavior in the T-dependence of log(r a ) as one approaches the glass transition, is 
independent of density and is thus an inherent property of each glassformer. This leads us to 
suggest a modification of the widely used Angell plot. 

As a more extensive data base on the P and T dependences of the a-relaxation time (and of 
the viscosity) of glassforming polymers and liquids now becomes available [3] , and reasonably 
accurate equations of state exist to convert P, T results to p, T results, the question arises of 
how to best organize the data, in a model-free and physically meaningful manner. In the case 
of molecular liquids, it has been recently suggested that the p dependence could be described 
through a single parameter, an effective activation energy E oc (p) characteristic of the high-T 
Arrhenius regime, i.e., log(r a (p,T)) = F(E oa (p)/T). The same procedure, however, is not 
applicable to glassforming polymers because most of them cannot be studied at T high enough 
for reaching an Arrhenius regime, so that E OQ (p) cannot be directly determined. In the cases 
studied in [4] E oa (p) was a quite featureless, increasing function of p (in the accessible range); 
rather than using an (uncontrolled) extrapolation of the data at high T, it then seems more 
sensible to choose a simple predefined function of p, e(p), with as few adjustable parameters 
as possible, to build a scaling variable X = e(p)/T. An appealing choice is a one-parameter 
power law, p x , x being material-specific. Such a choice has been used to represent a-relaxation 
time data for the fragile glassforming liquid ortho-terphenyl (o-TP), with x = 4 [5]; it is also 
consistent with a recent study of the density dependence of the glass-transition line T g {p) of 
several polymers (up to few kbars), where it was found that T g (p)ap x , with x « 2.0 for atactic 
poly (propylene), x m 2.7 for poly(styrene), and x w 3.4 for poly(carbonate) [6]. We have thus 
analyzed the a-relaxation time T a (p,T), obtained from dielectric relaxation measurements 
made on PECH, PVME, PVAc, and PMMA over a wide range of T and P (all experimental 
details are given in [7-9]) and converted to (p, T) via the equation of state proposed by Sanchez 
and Cho [10], as a function of p x /T, x being an adjustable parameter. As illustrated in Fig. 
1 , a very good scaling of the data is indeed reached. We find x w 2.7 for PECH and PVME, 
x ~ 1.4 for PVAc, and x ~ 1.25 for PMMA. We have applied the same analysis to the a- 
rclaxation time data obtained on 1,4-PB [11] with x « 1.8. In all cases a good data collapse 
is found with Eq. (1) and e(p) — p x /T : see inset of Fig. 4. 

Before discussing some of the consequences of the scaling, it is worth stressing that despite 
its appealing form, the power law behavior of e(p) may not convey much physical content. 
For o-TP it was suggested that the p 4 /T scaling was reminiscent of the well-known scaling 
property of simple liquids interacting via a soft repulsive pair potential, v(r) = e(^) n , with 
n = 12 as in the Lennard- Jones model [12]: indeed, for such systems the long-time properties 
depend on p and T only through a single combination pT _3// ", hence for n = 12, pT -1 / 4 , or 
equivalently, p 4 /T. However, such an interpretation implies a strong corollary, that the excess 
thermodynamic properties of the system (measured with respect to the ideal gas contribution) 
should also depend on the scaling variable, pT~ 3 ^ n or p n / 3 /T, only. We have checked that 
this is not the case for any of the polymers considered here (nor for o-TP): as illustrated in 
Fig. 2 for PVME and PB, the quantity (-^) obtained from the existing equations of state 
is not a unique function of p x /T with x determined from the scaling plot of the a-relaxation 
time. Not unexpectedly, in the available range of p, polymers do not behave as soft repulsive 
spheres. The power law p x should thus be considered as a convenient, one-parameter way 
of representing the p dependence of the activation energy scale oc e(p), a scale that, when 
data are available at high enough T, can be identified with the empirical Arrhenius activation 
energy E OQ (p) [4]. Actually, a good agreement to the a-relaxation time data can also be 
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Fig. 1 - Scaled plot of log(r a (p,T)) versus X = p x /T for PVME(i = 2.7) and PVAc(a; = 1.4). Data 
were obtained at atmospheric pressure and along different isotherms and isobars up to 300MPa. 



obtained by choosing a one-parameter linear description, e(p) = p — p* : see Fig. 3. Again, 
no special emphasis should be put on the fact that e(p) is a linear function of p; this merely 
reflects the fact that over the limited range of density that is accessible experimentally, e(p) 
increases monotonically with p in a featureless fashion. 

Besides the fact that it helps organizing all experimental data obtained by changing P 
and/or T in a simple and rational way, the main consequence of the scaling expressed by 
Eq. is that "fragility", which characterizes the super- Arrhenius T-dependence of the r Q 
and the viscosity, is an intrinsic property of a glassformer, in that it does not vary with 
density. The concept of fragility has proved to be most useful in the study of glassforming 
liquids and polymers and is now part of the very lexicon of such studies. Quantifying the 



4 



EUROPHYSICS LETTERS 



H 
Cl 



1.0 



0.8 



0.6 



0.4 



0.2 - 



0.0 



_3 

%^p=1.09g.cm 1Q 



0.6 (1< 

0.4 

0.2 



-3 PB 0.0 



p=1.02g.cm" 
0.04 0.05 0.06 

pi" - 55 



□ □ 



□ □ 



□ □ 



□ □ □ 



□ □ 



□ □ ° 



□ □□□□□n nD 



□ □□□□□□□□□□ 



PVME 



0.115 0.120 0.125 0.130 0.135 0.140 

pT -0.37 



Fig. 2 - Equation of state Pj pT versus X — p x /T for PVME calculated from [10] along isobars; 
inset : 1,4-PB along two isochores. Note the absence of scaling. 



degree of fragility of a system can be conveniently done without having recourse to fitting 
formulae by considering the steepness index at constant density, m p . T , defined as 
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Fig. 3 - Scaled plot log(r a (p,T)) versus X = e(p)/T with e(p) = p - p* for PVME (p* = 0.62) and 
PVAc (p* = 0.29); same symbols as in Fig.l. 



dlog(T a ) 

m P,r = .g. (T = T T ), (2) 

where T r is the temperature at which r a = r (expressed, say, in sec.) at the given p. By 
using Eq. 0) and defining A" r as the value of the scaling variable such that F(X T ) = log(r), 
it is easy to show that m p . T = X T F' (X T ), where F'(X T ) — dF Jx^ > 1S independent of density. 
Thus fragility, when properly defined, is an intrinsic property of the glassformer. Breakdown 
of this feature may occur when major structural changes take place in the system under 
compression, as for instance in tetrahedrally bonded liquids such as SiC>2 or H2O in which 
the local coordination is known to change with density. The empirical observation that fragility 
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Fig. 4 - Modified Angell plot : log(T a ) versus X/X g where X g is the glass transition value (r Q = 
lOOsec or r] — lQ ri mPa.s) for all polymers shown in inset. For comparison we also show as full 
and dashed lines the scaled curves for o-terphenyl and glycerol, respectively [4]. Inset : scaled plot 
log(r a (p,T) versus X = p x /T for several polymers, PVME, PECH(:r = 2.7), PMMA(i = 1.25), 
PVAc(:r = 1.4), 1,4-PB (as = 1.8) and o-terphenyl (x = 4). 



may vary with pressure is a consequence of the possible change in the relative contributions 
to slowing down due to T and to p. Indeed, the steepness index at constant pressure, TOp, T , 
denned by an equation similar to Eq. © is related to that at constant density through 

m P>T (P) = m P)T (l + a P /\<Xr\), (3) 

where ap = —p~ 1 (dp/dT)p and a T = —p~ 1 (dp/dT) T are the isobaric and isochronic 
coefficients of expansivity; the ratio ctp/\a T \ characterizes the relative effect of p over that 
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of T at constant P when r Q = r [13]. (1 + ap/\a T \) is also equal to [14], with Hp = 
{d\n{T a )/d(l/T)) P a,ndE v = (dln(r a )/d(l/T)) p ). One can see from Eq. ^ that the fragility 
measure at constant P is always larger than the constant-p fragility and that its variation with 
P depends on that of the ratio ap/\a T \] the amplitude and the sign of this latter is species- 
specific : for instance it has been shown to decrease with increasing P in the case of glycerol 
(resulting in a decrease of the fragility with pressure [15]) and to increase with P in the case of 
1,4-polybutadiene [13]. Combining Eq. and Eq. (j3J) gives ctp/\a T \ — ttp^wT T , which 
in the case of a power law behavior e(p) = p x simply reduces to ap.x.T T , where both ap and 
T T depend on P; usually ap decreases with P [16] and T T increases with P, lead to a variety 
of behavior for ap/\a T \ (see above). 

The fact that fragility is an intrinsic property of a glassformer is best represented by 
modifying the now standard Angell plot [2] in which log(r a ) is shown versus the inverse 
scaled temperature T g /T at constant (usually atmospheric) pressure : we suggest instead, 
when enough data is available, to plot log(r a ) versus X/X g , where X = e(p)/T is the scaling 
variable introduced above and X g its value when r a reaches a characteristic " glass transition" 
value, say r Q = lOOsec for dielectric relaxation data. The steepness of the log(r a )-vs-X/X g 
curve is a measure of the intrinsic fragility of a system, since indeed, according to Eq. QJ, 



dlog(r a 

m x - 



d(X/X g 



= X g F>(X g ), (4) 

x=x a 



is equal to m PiT at the chosen glass transition point. The modified Angell plot is shown in 
Fig. 4 where all the polymer data considered in this study (see inset of Fig. 4) are displayed 
together with results for a "fragile" liquid o-terphenyl and an "intermediate" one, glycerol [4]. 

In conclusion, we have shown that the P and T dependences of the a-relaxation time of 
glassforming polymers and liquids can be combined, after conversion to p—T data, in a function 
of a single scaling variable X = e(p)/T, where e{p) is well fitted by a power law p x , x being 
species-specific. On the practical side, this provides a simple way to organize and display the 
increasing amount of isobaric and isothermal relaxation time data on glassforming polymers 
and liquids. On the theoretical side, this suggests that fragility is an intrinsic property of a 
glassformer, which is best illustrated on a modified Angell plot log(r a ) versus X/X g , and that 
e(p) characterizes the density dependence of the "bare" activation energy scale. 
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